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Q\ ■ Abstract. The present article discusses the connection between 

, exactly-solvable Schrodinger equations and the Liouville transfor- 

mation. This transformation yields a large class of exactly-solvable 
potentials, including the exactly-solvable potentials introduced by 
Natanzon. As well, this class is shown to contain two new families 
of exactly solvable potentials. 



O 
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> ■ 1. Introduction 

The study of exactly solvable Schrodinger equations dates back to the 
\q \ very beginnings of quantum mechanics. As examples one can site the 

harmonic oscillator, the Coulomb, the Morse ||Mor29| , Poschll- Teller 



PT33], Eckart |Eck30|, and the Manning-Rosen |[MR33|1 potentials. 



One can argue that in each of these cases the exact-solvability comes 
about because the Schrodinger equations in question can be trans- 
formed by a gauge transformation and by a change of variables into 
either the hypergeometric or the confluent hypergeometric equation. 
^ ! To be more precise, in each of the above cases there exists a gauge fac- 

tor a(z; E), which depends on the energy parameter E, and a change 
of variables z = z(r), which does not, such that solutions to the corre- 
sponding Schrodinger equation, 

-ip"(r; E) + U(r)ip(r; E) = Eijj(r; E), (1) 

are of the form 

V>(r; E) = exp[a(z(r); E)\ 0(z(r); E), (2) 

where (f>(z;E) is either F(a, (3; 7; z), the Gauss hypergeometric func- 
tion, or $(«;7;^), the confluent hypergeometric function, and where 
the parameters a, (3, 7 are themselves functions of E. The just men- 
tioned types of special functions are well understood, and as a con- 
sequence one can explicitly calculate the bound state and scattering 
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information for the corresponding potentials. In light of these remarks 
the following question is of interest. 

Problem 1. Given a collection of functions, T = {<ft(z)}, find all pos- 
sible potentials, U(r), such that there exist an E -dependent gauge fac- 
tor, a(z; E), and an E -independent change of variables, z(r), such that 
the solutions of equation ([!]) are of the form shown in (pj) . 

For the cases of hypergeometric and confluent hypergeometric func- 
tions, Problem |I| was solved by Natanzon in ||Nat71|| . The correspond- 



ing classes of exactly- solvable potentials have come to be known as 
Natanzon's hypergeometric and confluent hypergeometric potentials, 
and have been the subject of some discussion in the literature ||Gin84| 



US93|| ||VVAG89|| . The purpose of the present article is to review Natan- 



zon's approach and then to enlarge Natanzon's class of exactly-solvable 
potentials by allowing to come from a larger class of special functions, 
namely the solutions of the following class of differential equations: 

A(z)<p"(z) + B{z)(j)'{z) + C(j){z) = 0, (3) 

where A (z) is a non-zero real polynomial of degree 2 or less, B(z) is a 
real polynomial of degree 1 or less, and C is a real constant. 

Prior to Natanzon, Problem |l| was considered by Bose in [|Bos64 



and by several other authors ||Man38|| ||BS64|| . Bose's paper is notewor 



thy because it introduced the approach that was followed by Natanzon 
in his classification. This approach relies on two techniques: a cer- 
tain canonical form for linear, second-order differential operators, and 
the Liouville transformation. The Liouville transformation will be de- 
scribed in Section ^| and the Bose-Natanzon approach in Section |3|. The 
solution of Problem [I] for the case where T is the set of solutions of 
equation (|3|) is given in Section f|. The resulting collection of potentials 
includes Natanzon's hypergeometric and confluent hypergeometric po- 
tentials, as well as two new classes of exactly-solvable potentials. These 
new potentials will be discussed in Section ||. 

2. The Liouville Transformation 
Consider a linear, second-order differential equation 

a(z)(f>"(z) + b{z)(f>'{z) + c{z)(f>{z) = 0, (4) 
Dividing through by a(z) and making the gauge transformation 

fa) = exp H' ^dt\ (f>(z) (5) 
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changes the equation into the following self-adjoint, canonical form 

4>"(z) + I(z)4>(z) = 0, (6) 
where the potential term is given by 

I = — - (4ac - lab' + 2ba -b 2 ). (7) 

Clearly, I(z) is an invariant of equation @ with respect to gauge trans- 
formations and multiplication by functions, and this is why equation 
(|6|) is being called a canonical form. Henceforth, I(z) will be called the 
Bose invariant of equation (f|). 

A change of the independent variable, say z = z(r), will transform 
equation (H) into 



[z'(r)]-Y(r) 



: ' (l ^-4>'(r)+I(z(r))t(r) 



0. 



where 



where 



[z'(r)} 3 

[z(r)). The corresponding canonical equation is 

i[)"(r) + J(r)i/)(r) = 0. 

ip(r) = [z'(r)]~2(f)(r), 
J(r) = [z'(r)} 2 I(z(r))+ 1 -{z J r} J 



(9) 
(10) 
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and where the curly brackets term denotes the Schwarzian derivative 
of z with respect to r, namely 

{z,r} 

The above process of going from one self-adjoint equation to an- 
other by means of a change of variables has been named the Liouville 
transformation in | Qlv74 1 , and the Liouville-Green transformation in 
|Zwi92|| . The Liouville transformation arises naturally in the context 
of WKB approximation (see Chapter 6 of ||Ulv74j| ) , and also underlies 
the following classical theorem due to Schwarz (see ||Hil76| 1, Theorem 
10.1.1 or ||(Jlv95|| , Theorem 6.28 .) 

Theorem 2. The general solution to the Schwarzian equation 

{z,r} = 2J(r) 

has the form z(r) = ip2{ r ) / ^i^) , where ^{f) andipi(r) are two linearly 
independent, but otherwise arbitrary solutions of equation (|8|) . 
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In particular, this theorem implies that for every potential J(r), there 
exists a change of variables z(r) such that the corresponding Liouville 
transformation takes the equation <fi"{z) = to equation (§). Therefore, 
one can relate any two equations of form @ by a Liouville transfor- 
mation. It is for this reason that Problem |1] must be formulated with 
the condition that z(r) not depend on the energy parameter. Without 
this restriction the problem would be uninteresting; one would get a 
criterion that would be satisfied by all possible potentials. 



3. The Bose-Natanzon approach 

The approach in question rests on the following reformulation of 
Problem [TJ 

Problem |l[A. Given a collection of functions, T = {<p(z)}, find all 
possible I\{z) > and Iq{z) such that for some E-dependent gauge 
factor, a(z; E) the solutions of 

4>"(z; E) + [h{z)E + I (z)) fa; E) = 0, (11) 

are of the form 

fa;E)=exp(a(z; E)) (j)(z; E). 

Indeed, suppose that Ii and Io satisfy the above set of requirements. 
Let z(r) be a solution of 

z'(r) = (h(z))^ , (12) 

From formula ( |T0l) it follows that a Liouville transformation of equation 
(|TT|) based on the change of variables z = z{r) yields an equation with 
potential term E — U(r) where 

UM _ h(z) -Ah{z)I';{z) + h{I[{z)f 

U[r) ~W) MhiW ' { } 

Furthermore, the corresponding eigenfunctions will have the form 

i/>(r;E) = (I 1 {z))* fa;E). 

Therefore U(r) satisfies the criterion imposed by Problem [l]. One can 
also reverse the above argument to show that given a U (r) demanded 
by Problem [l], one can produce an Ii and an I that satisfy the criterion 
of Problem |l]A. In other words, the two formulations are equivalent. 
The Bose invariant for the hypergeometric equation, 



z(l - z)4>"(z) + (7 - (1 + a + (3)) <j)'{z) - a(3 <f>(z) = 0, 
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is given by 

T(z) 
Az 2 {\ - z 
where 



'CO - ,,,„, , 2 



T(z) = (1 - (a - p) 2 )z 2 + (2 7 (a + /3 - 1) - Aap)z + 7 (2 - 7 ) 

Note that every polynomial, T(z), of degree 2 or less can be obtained 
from some choice of a, (3, 7 . Therefore, in order to solve Problem |T]A 
one must determine all possible I\{z) and Io{z) such that for all E there 
exists a T(z; E) of degree two or less in z such that 

It is clear that this condition is satisfied if and only if T(z; E) = 
R(z)E + S(z) where R(z) and S^z) are polynomials of degree 2 or less, 
and such that R(z) > in the domain of interest. The determining 
relation for z{r) follows from fll2j) ; it is 

z (r) 



Setting R(z) = r 2 z 2 + T\Z + r , calculating {z,r} and plugging the 
result into ( |T0| ) one obtains the formula for Natanzon's hypergeometric 
potentials: 

_ S(z) + l f r 1 -2(r 2 +r 1 )z 5 (r\ - 4r 2 r ) \ z 2 (l-z) 2 

V 4 i?(z) 7 i?(z) 2 1 ' 

These potentials describe the solution of Problem |l] for the case where 
T is the set of hypergeometric functions. 

It is well known that one can transform hypergeometric functions 
into confluent hypergeometric ones by a certain limit process. Natan- 
zon obtained his confluent hypergeometric potentials by applying this 
limit process to his hypergeometric potentials. The resulting family 
of potentials can also be considered as a solution of Problem 1, but 
the corresponding T is not the set of confluent hypergeometric func- 
tions, but rather the set of scaled confluent hypergeometric functions; 
namely cj){z) = $(a; 7 ; uz), where u is an extra scaling parameter. 
These functions satisfy the following scaled version of the confluent 
hypergeometric equation 

z(f)"(z) + ( 7 — ujz) <f)'(z) — ujacj){z) = 0. 

The corresponding Bose invariant is 

-uj 2 z 2 + 2uh - 2a)z + ^(2 - V) 
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where again every possible second degree polynomial can occur in the 
numerator as one varies a, 7, u. Hence, by the same reasoning as 
above, the criterion of Problem [l|A will be satisfied if and only if 

R(z)E + S(z) 



h(z)E + I (z) 



4z 2 



where R(z) and S(z) are polynomials of degree 2 or less, such that 
R(z) > in the domain of interest. The determining relation for z(r) 
follows from (|12"D ; it is 

z'(r) = -^=. 

VW) 

Setting R(z) = r 2 z 2 + r 1 z + r , calculating {z,r} and plugging the 
result into flIU|) one obtains the formula for Natanzon's confluent hy- 
pergeometric potentials: 



S(z) + 1 ( r\ 5 (r\ — 4r 2 ro) 



z 



2 



U = ^[zT + {-I-l R(z) - r2 )RW (15) 

4. Generalized Natanzon potentials 

The present section is devoted to the solution of Problem |l] for the 
case where T is the set of solutions, <f>(z), of equations of type (§). The 
Bose invariant for equation (H) is given by 

T( \ T(Z) 

I(z) 



A(z) 



where T(z) is a polynomial of degree 2 or less determined by quadratic 
combinations of the coefficients of A(z), B(z), and C. The exact for- 
mula for T(z) is not important; it can be readily recovered from equa- 
tion What is significant, is that for a fixed A(z), one can obtain 
every possible T(z) of degree 2 or less, from some choice of B(z) and 
C . Consequently, using the reformulation given by Problem [TJA, one 
must seek all possible I\{z) and Iq{z) such that for all values of E, 
there exist T(z; E) and A(z; E) both of degree 2 or less, such that 

/ lW £ + /„W = ^. (16) 

Lemma 3. Suppose that for all E there exist T(z; E) and A(z; E) such 
that relation ([II]) holds. Then, there exist E -independent polynomials 
A(z), R{z), S(z) of degree two or less such that I\ = R/A 2 and Iq = 
S/A 2 . 
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Proof. Setting E = in (|16D, one infers that I = Tq/Aq where the 
degrees of To and Aq are two or less. Similarly by setting E = 1 one 
infers that Ji = Tq/Aq — T\jA\ where the degrees of T\ and A\ are two 
or less. Thus, relation (|TJ|) may be rewritten as 

ET\ (1-E)T ET\A\ + (1 — E)TqA\ T(z;E) 
A\ + A\ A\A\ A{z-Ef [ ' 

Given rational functions Pq/Qq and P\/Q\ where the numerators and 
denominators are relatively prime polynomials, it's easy to show that 
the denominator of the reduced form of 

Qo Qi 

is the least common multiple of Qo an d Qi for all but a finite number 
of A values. This observation implies that the least common multiple 
of Ai and A must have degree less than or equal to 2. 

The rest of the proof will be done by cases, based on the degree of 
A± and Aq. If both Aq and A\ are constants then there is nothing to 
prove. 

Next, consider the case where one of Aq or A\ is a constant, but 
the other one isn't. Without loss of generality suppose it is A\ that 
is constant. Then T\ must be constant also, for otherwise the linear 
combinations of TqA\ and T\A^ would not always yield polynomials of 
degree 2 or less. Therefore the lemma is true for this case also. 

Suppose next that both Ai and Aq have degree 1. If A x and A 
have the same root, then there is nothing to prove. If A\ and A 
have different roots, then both 7\ and T must be constants, because 
otherwise linear combinations of TqA\ and T%Aq would not always yield 
polynomials of degree 2 or less. Hence one can write 

r T x Al _ TqA\ 

A 2 A 2 ' A 2 A 2 ' 

and this proves the lemma for the case under consideration. 

Finally, suppose that one or both of Ai and Aq is second degree. 
Without loss of generality assume that A\ has degree 2. Consequently, 
A\ must be the least common multiple of A and Ai, i.e. A must be 
a factor of A%. Now Aq cannot be a constant, because generically, the 
degree of linear combinations of T\A^ and TqA\ would be greater than 
2. If the degree of Aq is 2, then there is nothing to prove. The last 
possibility is that A\ = AAq, where both Aq and A have degree 1. In 
this case 

ETi + (1 - E)A 2 T 

hh + 1 = , 
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and hence T must be a constant. But one can therefore write 

T n A 2 



A\ ' 

i.e. the lemma is also true for this last case. □ 

Using the above Lemma, as well as the formulas ( |12D and (|13D one 
can now give the solution of Problem [l] for the case where JF is the set 
of solutions to equation (^). The desired potentials have the form 

_ -S(z)+V(A)/4 ( 3R"(z) 5 PGR) R'(z)A'(z) \ A{zf 
[ > R(z) ^\ 2 4i?(z) + A(z) J 4R(z) 2 ' { j 

where A(z), R(z), S(z) are polynomials of degree two or less, where 
V denotes the discriminant operator, and where z{r) is a solution of 

z'(r) = -4^=- (19) 

Henceforth this class of potentials will be referred to as the generalized 
Natanzon potentials. 



5. New exactly-solvable potentials 

Note that the presentation of the generalized Natanzon potentials 
given by equations (|T8| ) and ( Jl9| ) is invariant under affine substitutions, 
z i — ► dz ~\~ b. Consequently, no generality will be lost if one restricts 
A(z) to one of the following 5 possibilities: 1, z, z 2 , z(z — 1), z 2 + 1. 
The second and the fourth case yield, respectively, the Natanzon con- 
fluent hypergeometric, and the Natanzon hypergeometric potentials. 
If A(z) = z 2 , then the substitution z \— > 1/z transforms ( |TBD and (|TP|) 
into the corresponding forms for the case of A(z) = z. Therefore, if 
A(z) = z 2 , one again obtains the Natanzon confluent hypergeomet- 
ric potentials. With a bit of work one can check that at the level 
of solutions to the respective forms of equation (|3|), this transforma- 
tion corresponds to the to the well known identity (see chapter 6.6 of 
jEE53l ): 

2 F (a, a + 1 - 7; 1/z) = z a ^(a, 7; z). 

where \1/ is the confluent hypergeometric function of the second kind. 

The two remaining cases, A{£) = 1 and A(z) = 1 + z 2 , yield new 
families of exactly solvable potentials. These will be referred to, re- 
spectively, as case 1, and case 5 potentials, and will now be examined 
in some detail. In each case it will be convenient to rewrite equation 
(0) using a certain choice of adapted parameters. These adapted pa- 
rameters will be denoted by Greek letters, and the resulting equation 
will be referred to as the primary equation. 
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The solutions of the primary equations can be given in terms of hy- 
pergeometric functions, and have a natural dependence on the adapted 
parameters. The actual potential depends on a choice of R(z). The 
potential parameters — these will be denoted using lower case Latin 
letters — and the energy parameter, E, will turn out to be related to 
the adapted parameters by polynomial relations. It is important to 
note that for fixed potential parameters, and a fixed value of E one 
must solve these relations in order to obtain the corresponding values 
of the adapted parameters of the primary equation. 

An examination of formula flIB| ) shows that in order for U(r) to be 
non-singular, the ^-domain of the function shown in the right hand side 
of (|18| ) must not contain any roots of R(z). Singular potentials will 
not be discussed here, and this constraint greatly reduces the possible 
choices for R(z). 

According to equation (|19D the physical coordinate of the correspond- 
ing Natanzon potential is given by 



In all cases one can explicitly calculate the above anti-derivative, but 
the inverse, i.e. z(r), cannot in general be specified explicitly. Indeed, 
one may reasonably speculate that historically the study of Natanzon 
potentials was delayed by the fact that the inverse, z(r), of the above 
anti-derivative can be given in terms of elementary functions only for 
certain restricted choices of R{z) and A(z) (based on the remarks found 
in the first paragraph of |[Nat71| it would seem that Natanzon shares 
this viewpoint.) 

Nonetheless, a great deal of information about the inverse is avail- 
able. First, one can always take the domain of z(r) to be the whole 
real line; this is a consequence of the fact that one excluded those R(z) 
that have roots in the domain of z. One can calculate power series 
and asymptotic expansions for z(r) and use these as the basis for a nu- 
merical approximation. The graphs of potential curves that are given 
below were generated using this approach. 

In the subsequent discussion $ and ^ will denote the confluent hy- 
pergeometric functions of the first and second kind, and F will denote 
the usual hypergeometric function, 2 F\. For the source of this nota- 
tion, as well as the various properties of these functions the reader is 
referred to ||EB53|| . 
Case 1. A{z) = 1. 

Primary equation: 4>"{z) — 2u(z — (3)(j)'(z) — Aaui 4>{z) = 
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Primary solutions: 

0o = $(a; 1/2; u{z - (3) 2 ) , fa = (z-0) $(a + 1/2; 3/2; ^ - (3) 2 ) 

Note that the family of potentials under discussion is invariant under 
substitutions of the form z i— > -2 + k, k G R. In order to obtain a non- 
singular potential, R(z) must not have real roots, and consequently it 
is sufficient to consider the case R(z) = z 2 + 1. The most general form 
of the potential can then be obtained by a scaling transformation. 
Distance 1-form and physical variable: 



dr 



Vz 2 + 1 dz, r = ^ (^z^z 2 + 1 + \og{z + Vz 2 + 1)J 



z ~ ± V 7 ^, r -> ±oo. (20) 
n _az + b 3/4 5/4 



z 2 + l (z 2 + l) 2 (,2 2 + l) 3 
E = -u 2 , a = -2(3u 2 , b = u 2 {(3 2 - 1) -w(l-4a). (21) 

The resulting 2-parameter family of potentials are characterized by the 
presence of two wells separated by a barrier - - see Figure [|. The 
parameter a controls the degree of asymmetry; if a = 0, the potential 
is symmetric about r = 0. The parameter b controls the height of the 
central barrier. As b increases the wells become smaller; if in addition 
a = 0, then they disappear altogether. As b decreases the two wells 
merge into one. 
Eigenf unctions: 

fa = e -l/2-(-/3) 2 ( z 2 + • = 0; L (22) 

Bound states: 

A bound state occurs when —2a G N, and when uo > 0. This directly 
implies that there are infinitely many bound states if a ^ or if b < 0; 
and that there are no bound states, otherwise. 
Scattering. 

The scattering states occur when E > 0. Correspondingly, u> = iu, 
and a — ■? + ia, where Cj > and a are real. As a consequence ipo an d 
fa are real valued functions; this follows directly from the well known 
formula (see chapter 6.3 of [EB53] ): 

$(a; c; z) = e 2 $(c - a; c; —z). 

An asymptotically free eigenf unction, call it fa, can be given explicitly 
in terms of the confluent hypergeometric function of the second kind: 

fa = w « e -i/2^-/3) 2 ^2 + !/ 2; w ( 2 _ p f). 
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From the well-known asymptotic formula (see chapter 6.13 of | EB53|| ) 

\l/(a; c; z) ~ z~ a , z — ► +00, 
and from (^) it follows that 



as r — > ±00. Thus, asymptotically represents an almost free parti- 
cle traveling toward the center. The discontinuity in the direction of 
motion is caused by the fact that ^(z) is not regular at z — 0. The 
extra terms in the asymptotic phase appear because of the slow rate 
- on the order of r _1 an r" 1 / 2 , depending respectively on whether a 
is zero or not — at which the potential falls off toward zero. 

From the relation between confluent hypergeometric functions of the 
first and second kind (see chapter 6.7 of [|EB53|1 ) one obtains 



■00 = c ip{ + coVf) ^1 = sgn(z)(ciV>f + ciV>f). 

where 

cb = (-a;)-« , c 1 = (-c)4- F ® 



r(i-id)' v 7 r(|-id) 

It immediately follows that 

1 /Vo _ M = fr^ f , z -»• +00 

2 V c o ci y l^ f + i?^f 2^-00 
where the reflection and transmission coefficients are given by 

T = e ie (l - e- 2a7Ti )-\ R = e w {l - e 2ani )-\ 

where 

r(«) 

Case 5. A(z) = z 2 + 1 
Primary equation: 

(1 + ^"(z) + (i(l - 2p) + (2a + 1)^)0'^) + {a 2 - 5 2 )4> = 0. 

The above equation is related to the usual hypergeometric equation, 

C(i - C)0cc + (7 - C(« + + i))0c - = 0, 

by a linear change of parameters, and a complex-linear change of coor- 
dinates: 

a + j3 a — P a + (3 1 — iz 

o = , <5 = — — , p = 7 — , C = r, ■ 
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Primary solutions: 

01 = F (a + 5, a - 5; p + a; (1 - Lz)/2) , 

02 = i 71 + 5, a - 5; 1 - p + (x; (1 + iz)/2) 

To obtain non-singular potentials one must take R(z) without any 
real roots. The resulting family of potential depends of 4 parameters. 
A treatment of the most general potential, i.e. one depending all 4 
parameters would be unduly long, and not particularly illuminating. 
Thus, the focus here will on a more manageable 3 parameter subclass, 
namely the potentials that correspond to R(z) = z 2 + a 2 . 
Distance 1-form and physical variable: 



dr = — - — dz, r = sinh 1 ( - ) + y/l — a 2 tanh 



a 2 z 



z 2 + l W Vv/^T 

z ~ aexp ^Vl — a 2 tanh -1 — a 2 j j sinh(r), r — > ±oo. 



1 z 1 + 1 / z 2 + 1 \ 2 / 1 la? 



Potential: 

U ~ a 2 + z 2 + 4 z 1 + a 2 + {^z 2 + a 2 ) ^ 4 2 z 2 + 1 + 4 z 2 + a 2 ) 
E=-8 2 ,b = S 2 (l - a 2 ) - (p- - (a - , (24) 

These potentials fall off exponentially toward zero for large r. Setting 
p = \ one obtains potentials that coincide with a certain subclass of 
Natanzon hypergeometric potentials. The correspondence is given by 
the following substitution: z \— > 1 + z 2 , and at the level of solutions 
to the respective primary equations is described by the following qua- 
dratic transformation of the hypergeometric function (see chapter 2.11 
of fEBSp : 



F (a + S, a - 6; 1/2 + a; (1 - iz)/2) = F ((a + S)/2, {a - 5)/2; 1/2 + a; 1 + z 2 ) . 

The generic shape is that of two spikes for a > 1, and two wells for 
< a < 1; when a = 1 one recovers a modified Poschll- Teller potential. 
The parameter b controls the height/depth of the central spike/well, 
while c is the skew parameter that controls the degree of asymmetry 
in the potential. 

The case a > 1 results in the more interesting potential shapes, 
and thus will be the focus of the remaining discussion. Consider the 
symmetric potentials (c = 0) for a fixed value of a > 1. The number of 
extrema in the potential curve depends on the value of b. There are 3 
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critical values of b where the number of extrema changes: 

7a 2 - 3 + 3a- 2 -a 2 + 9 - 9a~ 2 , 3 

, , -a 2 + -. 

20 4 ' 4 

At these critical values of b some of the extrema merge, and one obtains 
some distinguished potential shapes; these shapes are shown in Figure 
0(a). 

At the value b = (—a 2 + 9 — 9a 2 )/4 the potential takes the form 



a 2 - 1 / ,z 4 ((a 2 -7),z 2 + 6a 4 -12a 2 
4a 4 V (z 2 + a 2 ) 3 



a 2 -7) 



One can show that z = r/a + 0(r 3 ) near r = 0, and hence the first 
3 r-derivatives of the potential vanish. As a consequence, one obtains 
a well with a very flat bottom. The potential also possesses two local 
maxima; these correspond to spike-like barriers on either side of the 
well. The variation of a in this type of potential shape is shown in 
Figure (b). The two barriers vanish precisely at the critical value of 
b = -a 2 + §. 

For c / one can obtain similarly distinguished potentials whenever 
b attains a critical value where the potential extrema merge. There is 
no exact formula for these critical values of 6; they must be solved for 
numerically. The resulting asymmetric potentials and their symmetric 
counterparts are shown in Figure |3|. 
Eigenf unctions: 



a 2 + z 2 \ 3 



V* = (1- iz) * + 2-4(i + iz) 2-2+4 {y^-t) 0i, i = 1, 2. (25) 
Bound states: 

An examination of relations (|24] ) will show that in order to obtain a 
potential with real coefficients, S must either be real or imaginary; and 
either a — | or p — | must be real, while the other must be imaginary. 
Note that the following transformation of the parameters is a symmetry 
of the potential: 

p^i(a-l/2) + l/2, an -i(p- 1/2) + 1/2. 

The transformation 5 — > —5 is also a potential symmetry. The presence 
of these two symmetries means that without loss of generality one can 
assume that a is real, that 9ft(p) = 1/2, and that 8 is either positive, or 
positive imaginary. With these assumption in place, ip2 is the complex 
conjugate of ipi, and the latter is square integrable if and only if 8 > 
and a + 5 G — N. After a bit of calculation one can show that this 
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criterion implies that for fixed a, 6, c, the bound states are indexed by 
natural numbers, N = — (a + 5) such that 

< + 2VWT^ - 1 

In particular, if c = (the symmetric potentials) then there will be 
no bound states if b > —1/4. If b < —1/4, then the number of bound 
states is equal to the largest integer smaller than 1/2 + v^— b. 
Scattering. 

The scattering states occur when E > 0, and hence without loss of 
generality 5 is positive imaginary. For reasons detailed above, a will be 
assumed to be real, while 9ft(p) will be assumed to be 1/2. To compute 
the reflection and transmission coefficients it will be useful to introduce 
two more solutions of the primary equation: 

03 = ( {iz - l)/2)- a - 5 F(a + 5,5 + 1 - p; 1 + 25; 2/(1- is)), 

4 = ((iz-l)/2)- (7 - 5 F(a-5,p-5;l-25; 2/(1 + is)). 



As per the formula in (p5|), let ip3 and ^4 denote the corresponding 
eigenfunctions. The usefulness of ^3 and ^4 is that they represent 
asymptotically free particles traveling, respectively, towards and away 
from the origin: 

^3 ~ K -5 e T(*+P+S-l/2) m/2 e -S\r\ ^ r _^ ±QQ ^ 

^4 ~ K'e^^- 6 - 1 ^ ^V |r| , r - ±00, 



where 



K = - exp ^ Vl — a 2 tanh 1 (\/l~ 



Relations between the regular eigenfunctions, ipi, 1P2, and the irregular 
ones, -03, -04 are given by (see chapter 2.9 of |[EB53|| ): 



^1 = C 3 ^3 + C 4 ^ 



'I- 



i'2 



— p ±7ri( CT +5)^ 3 _|_ _ e ±7ri( CT -5) 



c 4 e 



'4, 



where the ± in the second equation corresponds to the sign of z, and 
where 



c 3 



r(q + p)r(2(5) 
T(p + 5)T(a + 5) 



c 4 



i> + P )r(-2£) 
r( P -5)r(a-5)' 
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It follows that 



K 6 e -(a+8+p-l/2) Tri/2 / e 7Ti(a+5) _ ^2 




C 3 C 4 

r — > +00 

e dr + Re- 5r \ 
where elementary calculations will show that 

K^Yja - S)T{1 -a- 5)T(p - 5)T(1 - p - 5) 

2ixT(-25)T{l-25) 
Y /sin(7rcr) sin(7rp) i cos(7rcr) cos(7rp) \ 
\ sin(7r5) cos(7r5) / 
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Figure 1. Case 1: roles of the a and b parameters. 




Figure 2. Case 5 symmetric potentials: (a) critical val- 
ues of the b parameter, (b) variation of the a parameter 
in potentials with the middle critical b value. 




Figure 3. Case 5 asymmetric potentials: variation of 
the c parameter in potentials with (a) lowest critical b 
value (b) middle critical b value. 



